A GENERALISED RAYLEIGH- TAYLOR CONDITION FOR THE 

MUSKAT PROBLEM 

(^ . JOACHIM ESCHER, ANCA-VOICHITA MATIOC, AND BOGDAN-VASILE MATIOC 

^~>' Abstract. In this paper we consider the evolution of two fluid phases 

^ . in a porous medium. The fluids are separated from each other and 

also the wetting phase from air by interfaces which evolve in time. We 

reduce the problem to an abstract evolution equation. A generalised 

Rayleigh- Taylor condition characterizes the parabolicity regime of the 

problem and allows us to establish a general well-posedness result and to 

study stability properties of flat steady-states. When considering surface 

Q, ' tension effects at the interface between the fluids and if the more dense 

.^^ ' fluid lies above, we find bifurcating finger-shaped equilibria which are 

. , all unstable. 

-)— > 

1. Introduction 

The Muskat problem is a widely used model for the intrusion of water 
into oil sand. A linear analysis was performed in [21,22,24] where a relation, 
the so-called Rayleigh- Taylor condition, was found to determine two regimes 
\/^ , for the problem: a stable regime, when a flat interface is stable under small 

CN ' deviations, and an unstable one, when fingering occurs. 

If^ • Nonetheless, existence and uniqueness of classical solutions has been firstly 

^D I proven in [25] by using Newton's iteration method. In the last decade the 

problem has received more interest and was studied by means of complex 
analysis [23], energy estimates [2,4-6], power series expansions [16], or ab- 
stract parabolic theory [13]. These different approaches cover a wide spec- 
trum of questions related to the Muskat problem: local well-posedness, global 
existence of solution, singular solutions, stability properties of equilibria. 

It is worth noticing that all these papers mentioned above consider the 
situation when there is only one moving boundary, namely the one separating 
the fluids. Either one prescribes boundary conditions at two boundaries 
which are kept fixed during the flow or so-called far-field boundary condition 
are imposed. This setting corresponds to an abstract equation with only 
one unknown - the interface between the fluids. In the present paper we 
consider the more involved situation when there are two moving boundaries, 
one separating the two fluids and one separating the wetting phase from air 
(assumed to be at uniform pressure equal to zero). The fluids are located 
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in a porous medium (or a vertical Hele-Shaw cell) and are assume to fill 
together with the dry phase (air) the entire void medium. Moreover, we 
incorporate gravity and viscosity effects into the modeling as well as surface 
tension forces at both interfaces. The invertibility of a bounded operator 
permits us to re-write the problem as an abstract non-autonomous evolution 
equation 

dtZ = ^t,z), z{o) = Zo, 

where the variable Z parametrises both unknown interfaces. The tempo- 
ral variable t is induced into the problem by the boundary condition b for 
the pressure on the bottom of the cell. For this problem we find a gener- 
alised Rayleigh- Taylor condition in terms only of the boundary data 6, the 
viscosities fi± , and densities p± of the fluids of the following form 

bfi++gp+fi^ >0 and — !— — -{b - gp+) + g{p+ - p^) < 0, (1.1) 

/z+ -I- p^ 

which determines the parabolic character of the problem in the absence of 
surface tension effects. When including surface tension forces at both inter- 
faces we may drop condition (1.1). We steadily use in this paper the sub- 
script — for the fluid on the bottom of the cell and -|- for that above. After 
showing that the Frechet derivative dz^{0) generates a strongly continuous 
and analytic semigroup, parabolic theory provides local well-posedness of the 
problem and the principle of linearised stability may be applied to study the 
stability properties of the unique flat equilibrium which is determined for a 
fixed amount of fluid -|- (this quantity is preserved by the flow) and a certain 
constant boundary data. 

When considering surface tension effects at the interface between the fluids 
and the more dense fluid lies above we re-discover the global bifurcation 
branches obtained in [13, 14] which consist only of finger-shaped equilibria 
of the Muskat problem. The exchange of stability theorem due to Crandall 
and Rabinowitz [8] applies to this particular problem and we show that all 
small equilibria are unstable. 

The outline of the paper is as follows: we describe in Section 2 the math- 
ematical model and present the main results. Section 3 is dedicated to the 
proof of the well-posedness result Theorem 2.1, and in the subsequent section 
we analyse the stability properties of the unique flat equilibrium as stated 
in Theorem 2.5. In Section 5 we prove our third main result, Theorem 2.7. 
The calculations leading to the representation of ^^^(O) as a Fourier multi- 
plication operator are done in the Appendix. 

2. The mathematical model and the main results 

Let us start this section by presenting the mathematical model of the 
setting described in the introduction. Given ?ti G N and f3 G (0, 1) the small 
Holder space /i'""'"^(S) stands for the closure of the smooth functions C°°(S) 
in C"*"''^(S). We let S denote the unit circle and functions on S are identified 
with 27r-periodic functions on M. For later purposes we define h^ (S) as 
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Figure 1. The mathematical setting 

the subspace of /i"^+''(S) consisting only of even functions, h^ {^) is the 
subspace of /i'""'"P(§) consisting only of functions with integral mean zero, and 



K 



functions to be 



/in'r(S) := 



m+/3, 



S)n/i! 



e 



S). Furthermore, we define the set of admissible 



^:={/gC2(S) : I/I < 1/2}, 

Each pair (/, h) G ^^ determines two open and simply connected subsets of 
the porous medium, seen as S x ( — 1,2) C § x M, as follows: 

f](/):={(x,y) : -1< y < /(x)}, 
f](/,/i):={(x,y) : f{x)<y<l + h{x)]. 

Let r > and (/, h) : [0, T] — )■ U^ be given such that, at each time t G [0, T], 
the fluid — is located at ri(/(t)) and the fluid + at Q(f{t),h(t)) (see Figure 
1). The two fluids are assumed to be of Newtonian type and incompressible, 
and both interfaces are supposed to move along with the fluids. The problem 
is governed by the following system of partial differential equations: 



Am+ 


= 





in 


Q{f,h) 


An_ 


= 





in 


m), 


a,/,+ wm-9^^^ 


= 





on 


T{h), 


u+ 


= 


gp+{l + h) --jdKTih) 


on 


r{h), 


U- 


= 


b 


on 


r_i, 


U+ — U- 


= 


9iP+- P-)f + lwKr{f) 


on 


r(/), 


9tf+ ^^J d,u± 


= 





on 


r(/), 


/(O) 


= 


/o, 






/i(O) 


= 


ho 







(2.1) 



for t £ [0,T], where (/o,/io) ^ ^^ determines the initial domains occupied 
by the fluids. We used the variable / for parametrising the interface r(/) := 
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[y = /(^)] between the two fluids and T{h) := [y = 1 + h{x)] separates the 
fluid + from air. The unit normal v at r(/) [resp. T[h)\ is chosen such that, 
if T is the tangent, the orthonormal basis {r, z^} has positive orientation. We 
also write nvif) ^'^^ i^V[h) for the curvature of r(/) and T{h), respectively. 
Moreover 7^ [resp. 7^] is the surface tension coefficient of the interface 
separating the fluids from air [resp. the fluids]. 

The potentials u± incorporate both pressure and gravity force u± := p± + 
gp±y, with g the gravity constant. The velocity fields v±, which satisfy 
Darcy's law 

k ^ 

v± = Vu-t, 

^± 
are presupposed to be equal on the boundary separating the fiuid phases. 
Hereby, k stands for the permeability of the porous medium. On the fixed 
boundary r_i := S x { — 1} we prescribed the value of the velocity potential 
U-. For a precise deduction of (2.1) we refer to [11, 13,25]. 

Let a G (0,1) be fixed for the following. A pair {f,h,u^,u^) is called 
classical Holder solution of (2.1) if 

(/,/.) GC([o,r],v)nci([o,r],(/ii+"(s))2), 

u+(t) G 6uc2+"(0(/(t),/i(t))) and n_(t) G buc^^'' {n{f (t))) for t G [0,r], 

and if {f,h,u^,U-) satisfies the equations of (2.1) pointwise. We defined 
V := Vi X V2 to be the subset of U^ given by 

Vi := {/ G /j2+2sign(7™)+a(g) : f aU}, 
V2 := {h G /i2+2sign(7,)+a(g) ; h^U], 

where sign(O) = and sign(7) = 1 for 7 > 0. The space buc ''~"(il(/)) is de- 
fined as closure of the smooth functions with bounded and uniformly continu- 
ous derivatives BUC°^{n{f)) in BUC'^+'^{n{f)). The space buc^+°'{n{f,h)) 
is defined similarly. Moreover, since the potentials u± are determined, when 
knowing (/, h), as solutions of elliptic problems (see Section 3) we also refer 
to (/, h) to be the solution of (2.1). The first main result of this paper states: 

Theorem 2.1. Let {'^d:lw) G (0,oo)^ be given. 

There exist an open neighbourhood O of the zero function in (/i^"'"'^(S)) 
such that for all (/o,/io) £ O and b £ C([0, 00), /i2+"(S)) problem (2.1) pos- 
sesses a unique classical Holder solution T{-\ (/o, /iq)) defined on a maximal 
time interval [0, T(fo, /iq)) and which satisfies J-{[0, T{fQ, ho)); (/q, ho)) C O. 

The mapping {(t,/o,/io) : (/o, /^o) G 0,i G (0,T(/o, /iq))} ^ (/i^+"(§))' 

(t,/o,/io) ^J='{t;{fo,ho)) 
has the same regularity as b has. 

Remark 2.2. The conclusion of Theorem 2.1 remains valid if 7^ = or 
7^ = with the following modifications: if 7d = 7«; = we have to replace 
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(/i^+°(S))^ by (/i2+"(S))^ and require that b E C([0,oo),c + Oq) where Oq 
is a sniah neighbourhood of the zero function in /i^"*"" (§) and c G M satisfies 

CM+ +9P+fJ'~ > 0, (2.2) 

^^~^" (c-g/>+) + g(p+-p^)<0. (2.3) 

When (7d,7^) G {0} x (0,oo) [resp. (7^,7^) G (0, 00) x {0}] we replace 
(/i4+"(S))^ by /i4+°(S) X /i2+°(S) [resp. /i2+"(s) x /i^+a^g)] ^^^^ request that 
the constant c satisfies only equation (2.2) [resp. eq. (2.3)[. 

Relation (2.2) is a generalisation of the positive pressure condition imposed 
in [10-12[ to ensure well-posedness and stability of the one-phase Hele-Shaw 
problem without surface tension. Indeed, if the fluids have the same densities 
and viscosity, (2.2) re-writes c-|-5p+ > which is, up to a scaling, the same 
condition as in [10-12]. Moreover, it turns out that the Muskat problem 
without surface tension effects studied in [6, 13,26[ is similar to our problem 
if 7d > 0. Indeed, we have: 

Lemma 2.3. The volume of fluid + is preserved by the solutions of (2.1). 

Proof. The proof is similar to that of [10, Lemma 3.1]. D 

In order to establish similarity between our problem when 7d > and 
that in [13,26], we determine a special solution of (2.1) in the case when the 
volume of fluid -|- is equal to 2tt, i.e. 



fo-hodx = 0. (2.4) 

If initially /(O) = /i(0) = /o G K and b depends only on time, then 

r f(t) -L. 9P+-b 

fit) = --^ :. . ''- , /(0) = /o, (2.5) 



/^- f{t) + 



fj'++fj- 



and, by Lemma 3.1, f{t) = h{t) as long as the solution exists. If />+ = /9_ 
and ^_ > 11^ we obtain from (2.3) that if 7^ = 0, then b > gp+, thus /' 
is positive if /o is close to zero, meaning that the more viscous fluid drives 
upwards the less viscous one in the medium. This condition has been found 
also in [13, 26] to guarantee well-posedness of the Muskat problem studied 
therein. Moreover, if the Atwood number 

A,, := {fi+ -/i_)/(//+-F/i„) 

is zero, then (2.3) tells us that the more dense fluid must lay beneath in 
order to guarantee well-posedness of (2.1) when 7^ = 0, result similar to 
that in [6, 13[. 

Corresponding to the result in [13[, where an optimal value for the normal 
velocity at which water may replace oil in the absence of surface tension 
effects was found, we obtain herein an optimal value for the pressure on the 
bottom of the medium: 
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Remark 2.4. If the fluid below is water and that above oil, and we neglect the 
surface force at the interface between them, we find from (2.3) an optimal 
value 

Pmax ■■= g{p+ + p~) - g{p+ - P-)A~^ (2.6) 

for the pressure on the bottom of the porous medium below which water 
may drive upwards oil in a stable regime (no fingering occurs). 

Proof. Relation (2.6) is obtained form (2.3) in view of n_ = p- — gp- on 
r_i. The optimal value for the potential h is hmax = Pmax — 9P~, and if the 
boundary value b is close to this value we find that the solutions of (2.5) 
fulfill /' > 0, thus water drives oil upwards. This last assertion follows from 

gp+ - hmax = g{p+ - P~)AJ,^ < 

since it is well-known [3] that /j_|_ < /9_ and fj,^ > /.i_ (oil is less dense and 
more viscous than water). D 

We infer from (2.5) that if 6 = gp^ and /o = 0, then /(t) = h{t) = 
for all t > 0. Concerning the stability properties of the stationary solution 
(/, /i) = (0,0), which is the unique fiat stationary solution of (2.1) for b = 
gp+ and which satisfies (2.4), we state: 

Theorem 2.5. Let ^(i,'jw ^ [0,oo). Then: 

(i) If g{p-—p^)+'yw > 0, then the fiat equilibrium {f,h) = (0,0) of (2.1) 
is exponentially stable. More precisely, there exists positive constants 
M,6, andu such i/iaH/ ||(/o, /io)||/,2+2sign{7»)+c(§)xh2+2Bigi(7d)+a(g) < 6 
and {fo,ho) satisfies (2.4), then the solution {f,h) of (2.1) exists 
globally and 

||(/(t),/l(t))||^2 + 2sign(7™) + a(§)x/i2+2sign(7d) + <«(§) + 1 1 (^t /(*) ) ^^t ^(*) ) 1 1 (hl + " (§))2 
< Me "^ ||(/0,^0)||/i2 + 2Bign(7™) + a(§)x/i2 + 2sign(7d)+<i(§), Vf>0. 

(ii) If g{p- — p+) + 7uj < 0, then (/, h) = (0, 0) is unstable. 

Remark 2.6. When we study the stability of equilibria in Theorem 2.5 and 
Theorem 2.7 below we fixed a volume of fluid + equal to 27r, meaning that the 
initial data of (2.1) are presupposed to satisfy (2.4). This setting is imposed 
by Lemma 2.3, since the volume of fluid + is preserved by the solutions of 
(2.1). 

Theorem 2.5 is related to the exponential stability result established in [13, 
Theorem 5.3] for the Muskat problem with only one free boundary and is 
stronger than that in [16], where only stability is shown. Notice that if 
7u; = 0, then the flat solution is always stable, since p- > p+ is exactly the 
condition (2.3) which guarantees well-posedness of (2.1). Concerning the 
unstable case, numerical experiments [17] show that the interface between 
the fluids becomes very ramified, and dendrite like structures occur as time 
evolves if g{p- — p+) + 7^ < 0. 
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Figure 2. Steady-state solution of (2.1) when /)_ < /9+ and 
Iw >0. 



If 6 = gp+ and the volume of fluid + is equal to 27r, there exist also 
other stationary solutions of (2.1). They appear only in the unstable regime 
or sufficiently close to it, that is when 7«, > and the more dense fluid 
lies above in the cell. We show that for certain small 7io > there exist 
finger-shaped stationary solutions of (2.1), and therefore we shall refer also 
to (7^„, /, h) to be solution of (2.1). Given 1 < / G N, we define 

9{P+ - P~) 



ll ■-- 



P 



Theorem 2.7. Let b = gp+ and jw{p+ — p~) > 0. If (7^, /, h) is a stationary 
solution of (2.1) satisfying (2.4), then h = and (7«,,/) is a solution of the 
Laplace-Young equation 

7Mf) + g{p+-p^)f = 0. (2.7) 

The solution of (2.7) are, up to a translation, even and all even solutions of 
(2.7) can be represented as a disjoint union 

UfZi{Me)Ji{e)) : eGM}Uu£i(7m>7i)x{0}U(7i,oo)x{0}, 
with continuous functions 

(7/,/0:M^(0,oo)x{/E/i^+"(§) : ||/||c(§) < 1}, 1 < / G N, 

which, near e = 0, are real analytic and satisfy: 
^9{p+ - P^ " 



7«(e) =li + 
While 7i is even and 



-e^ + 0(e^), fi{e) = ecos{lx) + 0{e'). 



27r^g{p^ 



either \\fi{e)\\c[s) A^oo 1 or ||//(e)||c(§) A^oo 00. 

Additionally, the equilibrium (7^(5), /;(e),0) of problem (2.1) is unstable if 
\e\ is small. When I = 1 we have to assume e 7^ too. 
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Here B stands for Euler's beta function. Notice that the stationary solu- 
tions of (2.1), which satisfy (2.4) (see in Figure 2), are the same with the 
stationary solutions of the Muskat problem studied in [6, 13], where just one 
moving boundary is considered {h is chosen a priori to be zero). For a precise 
description of the global bifurcation branches (7;,/;) we refer to [14]. It is 
shown there that the situation ||/i(^)llc(S) /^e-s>oo 1 may occur only for small 
integers I. 

3. The evolution equation 

In order to solve problem (2.1) we re- write it as an abstract evolution 
equation on the unit circle. To do that we first transform system (2.1) into a 
system of equations on fixed domains by using the unknown functions (/, h). 
Let n_ := 0(0) and 0,^ := 0(0,0). Given (/, h) £ V we define the mappings 

4>f{x,y) ■■= {x,y + {l+y)f{x)), {x,y) G Q^, 

respectively (pf^h = {4>) ,h, 4>) ,h) ■ ^+ -^ ^{f,h) 

4'f,h{x, y) ■■= {x, y(l + h{x)) + (1 - y)f{x)), (x, y) G 0+. 

One can easily check that (j)f and (pf^h are diffeomorphisms for all (/, h) G V. 
These diffeomorphisms induce pull-back and push-forward operators (see e.g. 
[11]) which we use to transform the differential operators involved in system 
(2.1) into operators on the domains 0,± and their boundaries, respectively. 
Each pair (/, h) £ V induces linear elliptic operators 

A{f) : buc'^+'^{n.) -^ &uc"(0_), A{f)v := A(t>_ o c/yj^) o 0^, 
A{f,h) : buc^+^{Q+) ^ &wc"(fl+), A{f,h)v := A{v+ o 0"!) o ^f^^,, 

which depend, as bounded operators, analytically on / and h. Denote by 
tro the trace operator with respect to Fq := S x {0}. We associate problem 
(2.1) the following trace operators on Fq: 

B{f)v^ := kf,z' tro((V(T;_ o ct>/)\i-f', 1)) o <Pf), v. G 6^c2+°(0„), 

B{f,h)v+ ■.= kf,-Hro{{V{v+o<pjj^)\{-f',l))o^f^^), v^ G 6^c2+°(0+), 

which, seen as bounded operators into /i^~^"(S), depend analytically on / 
and h as well. Lastly, we define a boundary operator on Fi := S x {1}. 
Given {f,h) G V, we let 

Bi{f,h)v+ := -kf,-Hn{{V{v+o^j)^)\{-h' ,1)) o^f^f,), v+ G 6^c2+-(f7+), 

whereby tri is the trace operator with respect to Fi. 

With this notation one can easily verify that if (/, h, u+, m_) is a solution 
of (2.1), then (/, h,v+ := Uj^ocpf^h^V- := U-ocpj) solves the following system 



GENERALISED RAYLEIGH- TAYLOR CONDITION 



of equations: 












Aif,h)v+ 


= 





in 


0+, 




A{f)v. 


= 





in 


0_, 




dth 


= 


Bi{f,h)v+ 


in 


Ti, 




v+ 


= 


gp+{l + h) --fdK{h) 


on 


Ti, 




V- 


= 


b 


on 


r-i, 


* 


V+ — V- 


= 


9{p+- P-)f + lwi^U) 


on 


To, 




dtf + B{f)v. 


= 





on 


To, 




dtf + B{f,h)v+ 


= 





on 


To, 




/(O) 


= 


/o, 








hiO) 


= 


Hq 







(3.1) 



for all t G [0,T], where the transformed curvature operator k : /i^"'~"(S) — )■ 
/i2+"(S) is defined by k(/) := f" /{l + f^f/"^ . The notion of solution of (3.1) 
is defined analogously to that of (2.1). Notice that the parametrisation (/, h) 
is left invariant by the transformation above. In fact, one can see, cf. [11, 
Lemma 1.2] that each solution of (3.1) corresponds to a unique solution of 
(2.1). 

We introduce now solution operators corresponding to the system (3.1). 
Given / e Vi and (g,p) G /ii+"(S) x h'^+'^{^), we let T{f,q,p) G buc'^+°'{n^) 
denote the solution of the linear, elliptic mixed boundary value problem 

( A{f)v- = in Q_, 

B{f)v. = q on To, (3.2) 

V- = p on r_i. 

Further on, we define 5 : Vx(/i^+°(S))^ -^ 6itc ^"'"° ($7+) by writing 5(/, /i,p,r) 
for the unique solution of the problem 

A{f,h)v+ = in n+, 

v+ = p on Ti, (3.3) 

r;_l_ = r on Tq. 

It is convenient to write T{f,q,p) = Ti{f)q + 72if)p, where 

Tiif)q:={A{f),B{f),trr\0,q,0), T2{f)p := iA{f),B{f),tr)-HO,0,p), 

respectively S{f, h,p,r) = 5i (/, h)p + ^2 (/) h)r, with 

Siif,h)p: = {A{f,h),tr,tr)~\0,p,0), 

S2{f,h)r : = {A{f,h),tr,try\0,0,r). 

The operators Ti{f) and Si{f,h), i = 1,2, are bounded linear operators and 
they depend, in the norm topology, analytically on / and h too. 
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The key point of our analysis is the following observation. If (/, h, w+, w_) 
is a classical solution of (3.1) for to the initial data (/o,/io), then it must 
hold: 

(i) /(O) = /o and h{0) = ho, 

(ii) v.=Tif,-dtf,b), 

{in) v+ =S{f,h,gp+{l + h) - 'ydK{h),troV- + g{p+ - p-)f + -fwK{f)), 

(iv) dtf + B{f,h)v+ = 0, 

(v) dth = Biif,h)v+. 

Let us now show that from (ii) — (iv) we can determine the derivative dtf 
as a function of /, h, and t only. Indeed, we plug (ii) into (in) and (in) into 
(iv) to obtain the equation 

dtf + Bif, h)Si if, h) [gp+{l + h)- 7,k(/i)] 

+ B{t\ h)S2{f, h)[tTo T{f, -dtf, h) + g{p+ -p-)f + jMf)] = 0, 

which can be writen equivalently 

(id;,i+.(s) -Bif, h)S2if, h) tro TiUmf + B{f, h)S2if, h) tro T2if)b 

+ B{f, h)S{f, h, gp+{l + h)- jMh),9{p+ -P~)f + l^i^U)). 

(3.4) 
The linear operator which is evaluated at dtf is invertible, so that we obtain, 
by applying its inverse to (3.4), an equation expressing the derivative dtf in 
dependence of /, h, and t. Indeed, we have: 

Lemma 3.1. The set V contains an open neighbourhood W o/ with the 
property that 

G{f, h) := id^i+.(s) -B{f, h)S2{f, h) tro Ti{f) G £(/i^+"(S)) 

is an isomorphism for all (/, h) G W. 

Proof. The proof is based on a continuity argument. Namely, all the oper- 
ators defined in this section depend analytically on their variables and then 
so does Q too. Thus, it suffices to show that ^(0,0) is an isomorphism. To 
do that, we represent ^(0,0) as a Fourier multiplication operator. Given 
q € /i^+"(S) we let q = ^mq^^h{m)e^^^ denote its Fourier series expansion. 
A Fourier series ansatz yields for Ti{Q)q the following expression 



Ti{0)q{x,y) = t^{l+y)m + ^ E ^f^m^p-m, ^Me^ 



for (x,y) G Q^. Respectively, if r = X^^gg^C"^)^*"^^' then 52(0,0)r may be 
expanded as follows 

p2m p—my pTny 

52(0,0)r(x,y) = (l-2/)f(0)+ ^ ^2m _ ^ > Me'"^" 

meZ\{0} 
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for {x,y) £ Q^. Combining these two relations and taking the normal deriv- 
ative yields that 



G{o,o)q 

thus t/(0, 0) is an isomorphism. 



^H 



-q, yqeh'^^iS), 



D 



In virtue of Lemma 3.1, if the pair (/, h) maps into W, we may apply the 
inverse of G{f, h) to (3.4), and get 

dtf = ^i{tj,h), (3.5) 

with a nonlinear and nonlocal operator $1 defined by the relation 

cl>i(t, /, h) := - g-\f, h) {Bif, h)S2{f, h) tro r2(/)6 

+ B{f, h)S{f, h, gp+{l + h)- jMh),g{p+ -P^)f + iMf))- 

(3.6) 
Furthermore, from (ii) — (v) and (3.5) we obtain that h is solution of the 
equation 

dth = <^2{tj,h), (3.7) 

where the operator $2 is given by 
^2it,f,h):=Bi{f,h)S{f,h,gp+{l + h)-jMh),g{p+-p^)f + ln.K{f)) 

+ Biif, h)S2{f, h) tro r(/, -<l>i(i, /, h), h). 

(3.8) 
By Lemma 3.1 and relations (3.5)-(3.8) we found that all the solutions 
(/, h) of (3.1) which are contained in W solve the following abstract evolution 
equation 

dtZ = ^{t,Z), Z{^) = Zq, (3.9) 

where <l> := (<I>i, <^2) and we introduced the new variable Z := (/, h). Con- 
cerning the operator <I> we state: 

Theorem 3.2. The operator ^ has the same regularity as b, it is analytic in 
the variable Z, and if b{0) =: c G M, then df^i{0) and (9/i$j(0), i = 1,2, 

are Fourier multipliers with symbols {Xf{m))mez o-nd (A^(?Ti))m6Z) "^ = 1)2, 
respectively, given by 

k\m\ 



\{{m) := [A^(c - gp+) + g{p+ - P-) - 7«;"i^] 



(/x+ +/i_)tanh(|m|)' 



(3.10) 



\{{m) :-- 



Aj(m) := 



4p+ -P'-) + '^9P+P'- 
P+ + M- 

gp+p^ + cp+ 2 
: ^Idm 



9P- 



Iwm 



k\m\ 



{p^ + /i_)sinh(|m|) ' 
(3.11) 



fc|m| 



{p^ +^+)sinh(|m|)' 



(3.12) 
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At(m) := 



\ ^Idm 



k\m\ ii_ X^(m) 

' ' ^ 1^ ^ (3.13) 



/i4. tanh(|?n,|) /x_(. cosh(m) 



Proof. The regularity assertion is obvious. That the first order partial deriva- 
tives of $j,i = 1,2, with respect to / and h are Fourier multipliers follows 
from (6.6)-(6.9), relations proven in the Appendix. D 

We give now a proof of our first main result. 

Proof of Theorem 2.1 . We verify that the assumptions of [19, Theorem 8.4.1] 
are fulfilled by $. Theorem 2.1 is then a consequence of this result. For con- 
tinuity reasons it suffices in fact to show only that the derivative 9^<1>(0) 
generates a strongly continuous and analytic semigroup in C{{h}^l^{^))'^), 
i.e. 

-5z$(0) G ^(/l2+2sign(7™)+/3(g) ^ /t2+2sign(7d)+/3(g)^ (/,l+^(g))2) 

for some (3 G (0, a). By using the interpolation properties of the small Holder 
spaces 

{h"'' (S), h""^ (S))e = /j(i-''Vo+eai ^g)^ ^3^^4^ 

if ^ G (0, 1) and (1 — 9)aQ + Oai ^ N, we find then all assumptions of [19, 
Theorem 8.4.1] to be fulfilled. Here (•,•) denotes the interpolation functor 
introduced by Da Prato and Grisvard [9]. 

Let us first notice that derivative (9/j$i(0) maps /i2+2sign{7d)+/3^g-j contin- 
uously into /i^'^^(S) for some j3 G (0,/3). This property can be verified easily 
by using [11, Theorem 3.4], which is a multiplier theorem based on some 
generalized Marcinkiewicz conditions. Since by (3.14) 



/j2+2sign(7d)+/3(g) ^ (/il+/5(S), /jS+S sign(7d)+/3(g)) 



(1+2 sign(7d)+/3-/3)/(2 sign(7d)+l) 



we deduce, in virtue of Theorem 1.6.1 and relation (2.2.2) in [1], that c)^$(0) 
generates a strongly continuous and analytic semigroup exactly when 9j$i(0) 
and 9/i$2(0) are generators, i.e. 

- 9/$i(0) G ?^(/i2+2«'g^(7-)+/3(s),/ii+^(s))^ 

- 9^^2(0) G ?^(/l2+2sisn(7d)+/3(S)^/il+/3(S)). 

When considering surface tension effects this property holds independently 
of the boundary conditions, and when 7,,^, = and 7^ = this is true if 

\ {c-gp+)+g{p+- P~)<^ and ■ >0, 

respectively, with 6(0) sufficiently close to c G M in h'^'^"{'B). We refer to [11] 
where the generator property of a Fourier multiplier between space of peri- 
odic and continuous functions is explicitly verified when knowing its symbol. 
In virtue of (3.14) the proof is completed. D 
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4. Equilibria and stability properties 

As we mentioned earlier, if we consider a fixed volume of fluid + equal to 
2tt and if 6 = gp-\-, then (/, h) = (0,0) is the unique flat stationary solution 
of problem (2.1). Moreover, the reduced equation (3.9) is autonomous since 
^ does not depend on time for constant b. In order to study the stability 
properties of this equilibrium, as stated in Theorem 2.5, we shall use the 
principle of linearised stability, and need therefore to determine the spectrum 
of the derivative d^{0). Being a generator and taking into consideration 
that the small Holder space h'^~^°'{S) is compactly embedded into /i^^"(S), 
we obtain from [18, Theorem III. 8. 29] that its spectrum consists entirely of 
isolated eigenvalues with finite multiplicity. 

In virtue of Theorems 9.1.2 and 9.1.3 in [19] we know that the trivial 
solution (/, /i) = (0,0) is exponentially stable if the spectrum of 9$(0) is 
bounded away from the imaginary axis in the left half complex plane, and 
unstable if the infimum of the real part of all eigenvalues in the right half 
plane is positive. One can easily see that if A is an eigenvalue of d^{0), then 
it must be, for some ?7i G N, eigenvalue of the matrix 

X{{m) A^(m) 
A2(m) A2(m) 

where, for b = gp+, we obtained the simpler expressions for the multiplier 
symbols: 



X({m) = - [g{p^ - p+) + 7^m^] 
A^(m) = - [g{p^ - /)+) + Jn^m^] 
\'l{m) = - [gp+ +7d"^^^ 



km 



{p-^ + p^) tanh(?n) 

km 
{p-f- + P-) sinh(m) 

km 



{p- + /i+) sinh(?Ti) 



, ;, , , r nn km u„(cosh (TTi) — 1) + U-l cosh (m) 

A{^) = - [9P+ + Idm^] ^^^ ■ ', w s -^- 

^ ^ /i-l_ + p- p+ sinh(?n,j cosh^mj 

(4.1) 

Thus, the spectrum of d^{Q) consists only of the eigenvalues 



\{{m) + A5(m) ± ^/(A{(m) - A^(m))2 + A\'l{m)\{{m) 
A±(m) := ^- ^ , 

whereby m, G N. Easily, we see that A+(0) = 0, thus we find our selves in 
the critical case of stability when is an eigenvalue, which makes it difficult 
for us to establish the stability properties of the flat solutions. This is due 
to the fact that the volume of fluid + is preserved by the flow, and this 
property has not been included yet into our equations (3.9). We do this by 
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introducing a new variable f := f — h ^ h^ *''§°'^T'* "Z^' "(§). Then 

dtf = dtf - dth = <!>!(/, / - 7) - <I>2(/, /-/)=: ^2(/, 7), 

(4.2) 
9J = $i(/,/-/)=:^i(/,/), 

and with this new variable, problem (3.9) is equivalent to 

dtX = ^{X), X(0) = Xo, (4.3) 

where ^ := (^1,^2) and X := {f,f)- The trivial solution of (3.9) corre- 
sponds to the solution (/, /) = (0,0) of (4.3), so that we shall study the 
stability properties of the trivial solution of (4.3) which, as we shall see, it 
is more convenient. This since: 

Lemma 4.1. It holds that ^2(/>/) G ^o'''"(^) f^''" ^^^ if^f) ^'^ ^ -^^^ neigh- 

bourhood W C /i2+2sign(7. )+«(§) ^ ^2+2sign(7.-7»)+«(g)_ 

Proof. Given (/, /) as above, let again h = f — f. Setting p := gp+{l + h) — 
7dK(/i) and r := iiQT2{f)h + g{p+ — p-)f + jwK{f), we infer from (3.6) that 

cl>i(/, h) = -g^Hf, h)B{f, h)Sif, h,p, r), 

which can be reformulated as follows 

$i(/, h) = Bif, h)S2{f, h) tro TiimiU, h) - B{f, h)Sif, h,p, r). 

Using this relation and (3.8), we obtain 



^2{fJ)dx= / Bif,h)S{f,h,0,troTi{f)<^i{f,h))dx 
Js 

+ / Biif, h)S{f, h, 0, tro Tiimiif, h)) dx 



B{f, h)S{f, h,p, r)dx - [ Biif, h)S{f, h,p, r) dx. 
; Js 

Therefore, in order to prove our claim, it will do if we show that 

B{f,h)S{f,h,p,r)dx+ [Bi{f,h)S{f,h,p,r)dx = 

I Js 



for all (/, /i) G W and arbitrary (p, r) E (/i^"'""(S))^. Defining the harmonic 
function u^ := S{f, h,p, r) o <I>7^, we then have 



B{f,h)S{f,h,p,r)dx+ Bi{f,h)S{f,h,p,r)dx = / dnU+ds 

Js Js /^+ Jr(f) 

/ 9„ti+ ds = — / An_|_ d{x, y) = 0, 

JT{f,h) f^+ Jn{f,h) 

for the outward unit normal at drt{f, h), i.e. n 
and n = —v on r(/). D 



Xf) 
__k_ 

fJ'+ Jr{f,h) fJ'+ Jn{f,h) 

where n stands for the outward unit normal at drt{f, h), i.e. n = 1/ on T{f, h) 



5^(0) 
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We come now to the proof of our second main result: 

Proof of Theorem 2. 5. Studying the stabihty properties of the trivial solu- 
tion (/, h) = (0,0) of (2.1), under the constrain (2.4), is equivalent to study 
the stability of the trivial solution (/, /) = (0,0) of (4.3), where, by Lemma 
4.1, we have that 

^:>V^/ii+°(S) x/i^+"(S). 
From (4.1) we know that every component of the matrix operator 

" 5/^1 (0) 9j^i(0) " 

5/^2(0) 9j^2(0) 

is a multiplier with symbol (x({m))m.ez, (>'i{m))m&z\{o}, (^2("^))meZ, and 
(A2("^))m6z\{0} respectively, given by 

A{(m) = A{(m) + A^(m), X{{m) = -X'^im), 

A{(m) = A{(m) + X^{m) - A{(m) - A^(?n), A[(m) = -A^(m) + A^(m) 

for m / 0, A{(0) = —kgp-{ij^ + //_)~^, and A2(0) = 0. A simple computa- 
tion shows that the spectrum of the Frechet derivative d^ (0) , which coincides 
with its point spectrum from the same reason [18, Theorem III. 8. 29], is given 
by 

a(5^(0)) = {A±(m) : 1 < m G N} U {-kgp^{fi+ + fi^y^}. 
Let us start and estimate the eigenvalues A-|-(?tt,) with ttt. G N. Since [x i— t- 
2;/tanh(2;)] is increasing on [0, oo) we conclude that 

A_(m) < A{(1)/2 <0, Vm>l. (4.4) 

Let now g{p^ — p+) + 7^ > 0. We show that in this case also the other 
eigenvalues A+(?tt-) are of negative sign. Indeed, it holds that 

A f,^\ - _A A{(m)A^(m)-A^(m)A{(m) 

-(A{(m) + A^(m)) + y (A{(?ti) - A^(m))2 + 4A5'(m)A^(m) 

(4.5) 
and, in view of 

A{(m)A^(m) - A?(m)A|(m) = k^m^gp^ + ^,mp[g{p. -p^) + j^m'] ^ 

p+[p+ + P-) 

(4.6) 
we conclude that the spectrum of 9^(0) is bounded away from the negative 
half axis in M^. The assertion stated in Theorem 2.5 {i) follows at once, with 
the remark that the constant u found there may be chosen arbitrarily in the 
set (0,-max{A : A G cr(a^(0)}). 

On the other hand, if g{p- — p+) + ^w < 0, one can easily observe that 
A_|_(l) > 0. Moreover, since A±{m) — ?• —00 as m — ?• 00, we conclude Theorem 
2.5 (ii). This finishes the proof. D 
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5. Finger-shaped equilibria 

This last section is dedicated entirely to the proof of Theorem 2.7. If the 
tupel (/, h, u^,u-) is a stationary solution of (2.1), it must hold that u_ = b, 
U-^ is constant, and 

lwK{f)+g{p+- p-)f = u+-b (5.1) 

-ydK{h) - gp+h = b - u+. (5.2) 

on S. Equations (5.1) and (5.2), called Laplace- Young or capillarity equa- 
tions, have been studied intensively (see [15] and the literature therein) sub- 
jected to certain constrains at a fixed rigid boundary. Though, when dealing 
with periodic solutions, we easily get, cf. [12], that h must be constant also in 
the spatial variable, and if p^ < /j_ then also / is constant. Whence, equa- 
tions (5.1) and (5.2) may have nontrivial solutions {f,h) ^ E? only when 
7i„ > and p+ > p-. 

We are interested to determine only the steady-states (/, h) of (2.1) when 
b = gpjf- and which satisfy (2.4), i.e. Q(f,h) encloses the same volume of 
fluid as il+. Assume by contradiction that h = c for some c 7^ 0. Since 
6 — u_l_ = —gp^c, we get that 

iw^if) + g{p+ - p^)f = gp+c. 

On one hand, if / is constant, it must hold that / = p^c/{p-^ — p-), 
which contradicts (2.4) for c 7^ 0. On the other hand, the function p := 
f — p^c/{p^ — P-) solves the equation 

P" 

^"' (l+p/2)3/2 + 9iP+ -P-)P = 0. 

The solutions of this equation are, up to a translation, odd. Indeed, since 
p is periodic and nonconstant, it must hold p{xq) = for some xq € M. By 
translation, we may take xq = 0. The pair (p, q := p') is a global solution of 
the initial value problem 

_g(/>+-P^)^(l^g2)3/2 ' ( ? ) ^°) = ( p'(0) 

This is also true for {p,q){x) := {—p{—x),q{—x)), x € M. Whence, as we 
claimed, p is odd, so that Lfdx = p+c/(p+ — p-), which contradicts again 
c 7^ and (2.4). Consequently, if (jwjfjh) is a solution of (2.1) and (2.4), 
then h = and (7^1,/) solves the problem (2.7), which implies in turn that 
/ has integral mean equal to and an even translation by [14, Theorem 3.2]. 
Combining that particular result with Theorem 6.1 in [13] we obtain all the 
claims of Theorem 2.7 excepting the stability assertion. 

In the remaining part of this section we prove that the steady-state solu- 
tion (7i(e), //(e), 0) of problem (2.1) is unstable provided that e is sufficiently 
small. We rediscover first the global branches (7;,/z,0), 1 < I € N, at least 
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locally near (7^,0,0), by applying the theorem on bifurcations from simple 
eigenvalues, due to Crandall and Rabinowitz [7, Theorem 1.7], to the prob- 
lem 

where ^ is the mapping defined by (4.2). We shall refer to (7^^,, 0, 0), 7^„ > 0, 
as being a trivial solution of (4.3). In order to use 7^„ as a bifurcation 
argument we establish first analytic dependence of ^ on 7^^ and take the 
restriction 

* : (0,00) X We C ht+"{§) X /i^+^''^^(^'')+"(S) ^ hl+^{S) X /i^+"(§), (5.3) 

where We := W n (/if+"(S) x /io+^''*5''^^'*^+"(§)). That ^ is well-defined be- 
tween these spaces follows by using elliptic maximum principles and Lemma 
4.1. When 7«, ^ {7/ : 1 < ^ G N}, we infer from (4.4)-(4.6) that ah the 
eigenvalues of d,r r-,^(7u;, 0) are different from zero, thus d,r 7x^(7^;, 0) is an 
isomorphism. The implicit function theorem ensures that (7^„,0,0) is not 
a bifurcation point of the trivial solution. Otherwise, if 7^ = 7; for some 
I E N, then 



Ker 5 



,r 7'^^{^l,0) = span{(cos(/x),cos(/x))}. 



whereby (cos(/x), cos(/x)) is the eigenvector corresponding to the eigenvalue 
A_|_(l,7;) = 0. We considered that A± = A±(m,jy^), i.e. A± depends not only 
on m, but also on 7^. Also, the codimension of the image Im(?(., T^^(7;,0) 
is one since 

(00 00 \ 

Y+Yl Omcos(mx), Y^ bn,cos{mx) G /i^+"(§) X /lJ+"(S) 
m=l m=l / 

belongs to Inid,, r-.^i'^i, 0) if and only if 

/i+(cosh^(l) — cosh(Z)) -I- ^„(cosh^(Z) — 1' 



{ai,bi) = X{l 



/i+ cosh(Q 
for some A G M. Moreover, one can easily check that the mixed derivative 



d 



/ 1 _ cosh(Z) 

, „ r-,\I'(7u;, 0)(cos(/a;),cos(Zx)) = A I cos(/x), —7-; — cos(/x) 

'^J'J) \ cosh(/j 



whereby A = — 2/c/^7; ((/U+ -|- //_) sinh(/))^ , does not belong to the image of 
d,r r-,^(7j, 0). We conclude by [7, Theorem 1.7] the existence of a bifurcation 
curve 

(7^/^^) : {-Si,Si) ^ (0,00) X ht^-{8) X /i^+2^'^"(^'')+"(S) 

consisting only of stationary solutions of (4.3). Since they correspond all to 
a volume of fiuid -|- equal to 27r, it follows that 7; and /; = /; are, up to a 
parametrisation, restrictions of the functions obtained in Theorem 2.7. 

The stability properties of the equilibrium (7;(e), /;(e), 0) for (2.1) un- 
der the constrain (2.4), are equivalent with that of the steady-state solution 
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(7/(e), /j(e), /;(e)) of problem (4.3). For our purposes, Theorem 2.7, it suf- 
fices in fact to show that (7z(e), fi{£), fii^)) is an unstable stationary solution 
of the abstract Cauchy problem 

dtX = ^{-,^,X), X(0)=Xo, (5.4) 

where ^ is the restriction (5.3). 

Indeed, if e small and / > 2, then (7i(e), /^(e), fii^)) is an unstable solution 
of (5.4) since the eigenvalue A_|_(l,7;) of d,, 7;^(7;,0) is positive. 

For the stability of the stationary solution (71(e), /i(e), /i(e)), when \e\ is 
small and \e\ 7^ 0, it is important how the eigenvalue A_|_(l,7i) (which is equal 
to 0) perturbs for small e. Our main tool is the exchange of stability theorem 
[8, Theorem 1.16] due to Crandall and Rabinowitz. The assumptions of this 
theorem are satisfied by ^ since: 

(a) d,r r-|^(7i, 0) is a Fredholm operator of index with a one-dimension- 
al kernel; 

(b) a^^ (^j)^(7i,0)(cos(3;),cos(x)) ^ Im 9(^^)^(71,0); 

(c) (cos(x),cos(x)) ^ Im9.„ 7;^(7i,0). 
Letting J' denote the inclusion 

/i^"(S) X /i^J'^'s"(^'')+"(S) ^ /ii+°(S) X /iJ+"(S), 

in the terminology of [8], (a), {b), and (c) mean that is a 9 ,, ri^(7i,0)- 
simple eigenvalue and a ^-simple eigenvalue of d,r 7x^(71, 0). By choosing 
61 sufficiently small, we obtain from [8, Theorem 1.16] four continuously 
differentiable functions A : (7^ — 5i,7i + Si) — >■ M, /i : {—5i,5i) — )■ M, 
u:{-f,- 6i,-fi + 5i) ^ hi+»{S) X /i^+^^'S'^(^^)+"(S), and v : (-duSi) ^ 
/if+-(S) X /i^+'"S"(^'')+"(S) such that: 

9(^j)^(7«;,0)u(7^) = A(7^„)u(7^„) for 7^„ G (7^ - 61, j^ + Si), 

d(^fj)'^{liie),fi{£)Ji{e))w{e) = /u(e)u;(e) for e G (-(5i,(^i), 

A(7i) = /u(0) = 0, and ti(7i) = w{0) = (cos(x), cos(x)). 
Moreover, A'(7i) 7^ and 

Since A(7„,) is an eigenvalue of 0, , 7-, ^ (7^; , 0) and A(7i) = we get, by 
continuity, that A(7u,) = A+(l,7tt,) for all |7m,— 7i| < ^i- Moreover, A+(l,7u,) 
is positive for 7^^ < 7i, and negative if 7^ > 7i) thus A'(7i) < 0. In order 
to determine the sign of the eigenvalue /u(e), which is the perturbation of 
the eigenvalue of 5, , 7,^(71, 0), we need to specify the sign of 7x(e). From 
Theorem 2.7 we obtain in view of 7^(0) = and 7i(0) > 0, that 7i(e) and 
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£ have the same sign, thus /x(e) is a positive eigenvalue, and we are done 
by [19, Theorem 9.1.3]. 

6. Appendix 

We end this paper with a detailed proof of Theorem 3.2. Since the dif- 
feomorphisms used in Section 3 to transform the original problem (2.1) into 
(3.1) are given explicitly in terms of / and h, we obtain by direct compu- 
tation the following expressions for the elliptic and trace operators defined 
therein: 

i2 



^ _ 2^1 + yKJL + (il+illl , _J__\ l_ 

dx^ 1 + / dxdy \ (l + f)^ (l + ffjdy 

(1 + /)/"- 2/2 d 



(1 + y)- 



(l + /)2 dy' 






9x2 l + h-f dxdy {l + h-fY dy^ 

yh" + (1 - y)f" (h' - f')iyh' + (1 - y)f)\ d 



l + h-f {l + h-fY ) dy' 

^^^^~[TVT'''^y-^'''Tx 

^^f^^^~[TTh^'''d-y-^'''d^)^ 

Siif, h)= — — tri — - /i tri — 

11+ \\ + h — f oy ox 

In the following we take 6(0) = c G M and show that 5/$i(0) and dh^i{^), 
i = 1,2, are Fourier multiplication operators. We shall also determine the 
symbol of these operators which are the main ingredients when proving the 
Theorems 2.1 and 2.5. 

6.1. The derivative df^i{0). Let us start by determining the symbol of the 
operator 9/$i(0). Putting h = in (3.6) yields that 

$i(0, /, 0) := - g~Hf, 0)^(/, 0)52(/, 0) tro T2{f)c 

- G-\f, 0)B{f, 0)52(7, OMp+ - p.)f + j^Kif)] 

-g-i(/,0)^(/,0)5i(/,0)5P+. 

We show first that the derivative of J^ := Q~^ with respect to / in is a 
Fourier multiplier. Indeed, by the chain rule we have, due to J-{f, 0)Q{f, 0) = 
id;,i+.(S), that dfTiO,0)[f]M = -g-\0,0)dfg{0,0)[f]g-\0,0)M for all 
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/ G /i2+2sign(7„ )+«(-§) ^^^ M £R. So, we need to determine 

9/^^(0, 0)[/]M = - 5/^(0, 0)[/]52 (0,0) trori(0)M 

- B{0, 0)5/52(0, 0)[/] tro Ti{0)M (6.1) 

-e(0, 0)52(0, 0)troari(0)[/]M. 

Using the expansions found in the proof of Lemma 3.1, we find that 7i(0)M = 
fj.^k-^M{l+y) and52(0,0)trori(0)M = ^_A:-^M(1 - y), hence 

- 6(0,0)52(0,0) tro ari(0)[/]M = ^^/. (6.2) 

/^+ 

Concerning the first two terms of (6.1) we determine first an expansion 

for ari(0)[/]M and 5/52(0,0)[/]M. We start with 9/52(0, 0)[/]M. Elhptic 

estimates yield that the function 9/52(0, 0)[/]Af is the solution of Dirichlet 

problem 

Aw = -dfA{0,0)[f]S2{0,0)M = -M{l-y)f" in Q+, 

w = on Fi, 

w = on Fq, 

and, as in the proof of Lemma 3.1, we get that 

,_^ / pmy _ 2m -my \ ^ 

9/52(0,0)[/]M = - Y. ^( ^I^iT^i + {l-y)\f{m)e^^\ 

meA{o} 
Whence 

-6(0, 0)9/52(0, 0)[/] tro ri(0)M = M^ ^ ^IL^/(^)e-- - ^^/. 

(6.3) 
Consider now the function 97i(0)[/]M, which is the solution of 

Aw = -9^(0)[/]ri(0)Af = /i_fc-iM(l + y)/" in [7_, 

dyw = -fi-k^^dl3{0)[f]Ti{0)M = n_k^^Mf on Fq, 

w = on F_i. 

Expanding 

f=Yl /(»^)e'"" and w{x,y) = ^ u;^(y)e*'"^ 

mSZ mSZ 

we find that Wm is the solution of the following problem 



w 



," ^1 



m?Wjn = —Mfi^k ^f{7n)m^{l+y) in — 1 < y < 0, 



<(0) = Mfi_k-'f{m), 

Wm{-1) = 0, 

which has the solution Wm{y) = M fi^k^^ f {m){l + y) for all m G Z. Whence 
^Tm[f]M = M^i_k-\l + y)f, 
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and, with the convention that r7i/tanh(?7T,) = 1 if ?7i = 0, we determine that 

-5(0, 0)52(0, 0)tro5ri(0)[/]M = ^^^ ^ ^JL^/(^)e--. (5.4) 

/!+ ^-^ tanh(?7ij 

The relations (6.1)-(6.4) fuse, in view of Lemma 3.1, to 

^ ^ ' ^ (/i_+A^+)2 A.^tanh(m)-'^ ' ^ ' 

for aU / e /i2+2sign(7„)+a(g) g^^^j ^^ ^ ]^_ jj^ ^^.^gj. ^^ determine the deriv- 
ative (9j<I>i(0) two more steps must be done: we must find the expansions 
of the derivatives 5/72(0) and 0/5i(0,0). Since 72{f)M = M, we get that 
972(0) [/]M = 0, and proceeding similarly as we did before, yields 

,_^ / p-my _ p2m -my \ ^ 

dfSr{0,0)[f]M = Y. ^{ e^^-l + U " ?/)) /("^K'""- 

mez\{o} ^ ^ 

Combining all these relations, we finally obtain for df^i{0,0) that 

df^mif] = E A((m)/(m)e'-^ (6.6) 

mez 

whereby (A{(m))mgz is given by (3.10). 

6.2. The derivative 0/^2(0). The expansions found in the previous subsec- 
tion are very useful when studying 9/$2(0,0), since by (3.8) we have 

cl>2(0, /, 0) ■.=B,{f, 0)5i(/, 0)gp+ - Biif, 0)^2 (/, 0) tro ri(/)<I>i(0, /, 0) 
+ i3i(/,0)52(/,0)tror2(/)c 
+ Bi{f,0)S2{f,0)[j^^{f)+g{p+-p.)f]. 

In view of ^i(O) = k{c — gp^)/{fi^ + /t^-) we obtain the following represen- 
tation 

5/^2(0)[/] = Yl A{(m)/(m)e*™^ (6.7) 

meZ 

with symbol (A2(wi))meZ given by (3.11). 

6.3. The derivative dh^i{0)- In order to determine a representation of the 
Frechet derivative 5/,,<I>i(0) we have to investigate first the partial derivatives 
with respect to h of the solution operators defined in Section 3. By (3.6), 
we get for / = that 

$1(0, 0, h) := - g~\0, h)B{0, h)S2{0, h) tro T2(0)c 

- g-i(0, /i)e(0, /i)cSi(0, /i)bp+(l + h)- 7dK(/i)]. 
Differentiating the relation J-"(0, /i)^(0, h) = id/ji+c<(§), yields at /i = that 

a/,j-(o,o)[/i]M = -g-^(o,o)a;,g(o,o)[/i]g-i(o,o)M 
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for all M £ R. We infer from the definition of G that 

dhg{0,0)[h] = -dhBiO,0)[h]S2{0,0)troTi{0) - B{0,0)dhS2iO,0)[h]tToTiiO). 

Given h G /i2+2sign(7d)+a(g)^ ^Yie partial derivative 9/^52 (0,0)[/i]M is the 
solution of the linear Dirichlet problem 

Aw = -dhA{0,0)[h]S2{0,0)M = -Myh" in n+, 

w = on Fi, 

w = on Fq, 

thus dhS2{0,0)[h]M expands as follows 

dnS2{m[h]M= Y. ^n pm_p-^ -y]h{m)e-\ 
mez\{o} ^ ^ 

and similarly 

m6Z\{0} ^ ^ 

It follows then easily that 

dnHmm'i = . ^-^+^. , Y. -^^^(^y""""- 

(/i+ + /i_)^^smh(m) 

By definition, 5i(0,0)p = (A, tr, tr)-^(0,p,0) for all p G /i2+"(S), and one 
can easily check that if p = X^,„g2P(W')e*™'^, then 

,^^ p-my _ p-my 

5i(0,0)p = yp(0)+ Y em_p-^ v{m)e'^\ 

meZ\{0} 

Summarising, we obtain that 

dn^M\h\ = Y >^'l{rn)h{m)e"^- , (6.8) 

meZ 

with symbol (A5*(?7i))meZ given by relation (3.12). 

6.4. The derivative a/,$2(0,0). Putting / = in (3.8) yields 

$2(0, h, 0) =Bi(0, /i)5i(0, h)[gp+{l + h) - 7dK(/i)] 

- Bi(0, /i)52(0, /i) tro ri(0)$i(0, /i, 0) 

+ Si(0,/i)52(0,/i)tror2(0)c. 
Using the relations derived above, finally yields 

dh<^2mh] = Y A{m)h{m)e'^\ (6.9) 

meZ 

with symbol (A2(m))mgz given by relation (3.13). 
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